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^ \ Abstract: Gluinos in the mass range ~ 1^ — 3^ GeV are absolutely 

excluded. Lighter gluinos are allowed, except for certain ranges of lifetime. 
Only small parts of the mass-lifetime parameter space are excluded for larger 

^ ' masses unless the fifetime is shorter than ~ 2 x 10- fa) sec. Refined 

mass and lifetime estimates for i?-hadrons are given, present direct and indi- 
rect experimental constraints are reviewed, and experiments to find or defini- 
tively exclude these possibilities are suggested. 
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1 Introduction 



Short-lived gluinos could be defined to be those which decay before in- 
teracting hadronically in a detector or beam-dump. Their decay produces 
the lightest neutralino (lsp)| which escapes the detector or dump without 
interacting, carrying with it much of the gluino's energy and momentum. 
Short-lived gluinos are excluded for masses < 160 GeVflj] by the absence of 
characteristic missing energy events in the FNAL collider. Thus to be lighter 
than this, gluinos must be long-livedP]. It is natural for gluinos to be much 
lighter than squarks if their masses are entirely radiative in origin. In that 
case, if the SUSY and elecroweak symmetry breaking scales are less than 
~ 10 TeV, gluino and lsp masses will range from of order 100 MeV to of 
order 30 GeV[|], H. This is the mass range explored here. 

A gluino in the mass range ~ 1.5 — 3.5 GeV is excluded, whatever its 
lifetime, from the absence of a peak in the photon energy spectrum in 
radiative Upsilon decay. This is because two gluinos with mass in that 
range would form a pseudoscalar bound state, the r]g, whose branching frac- 
tion in T — > jr/g can be reliably computed using perturbative QCD and is 
predicted || |7|, || to be greater than the experimental upper bound |)|, p~0HP|. 

2 Generally, a superposition of the SUSY partners of the photon, Z° and neutral Higgses, 
not considering effects of a possible light gravitino. 

3 The long-lived gluino window was first pointed out in ref. [J2|. It was subsequently 
discussed in the work of ref. J3| . 

4 The range excluded by the CUSB experiment is incorrectly claimed to extend to lower 
gluino masses, by using the pQCD results of refs. |6|, J7|, || out of their range of validity. 
A detailed analysis of the actual excluded range in given in ref. |11| . The lower limit 
for validity of a pQCD, non-relativistic potential model description of an r\ g was taken to 
be ~ 3 GeV, mainly by analogy with the success of the same description of charmonium. 
However since the effective value of the coupling is so much stronger due to the larger 
color charge of the gluino in comparison to a quark, even a 3 GeV r]~ g may not be in the 
perturbative regime, in which case the range of validity of the CUSB procedure may not be 
even this large. Note that any gluino whose lifetime is longer than the strong interaction 
disintegration time of the rjg, i.e., r > ~ 10~ 22 sec, will produce the requisite bump in the 



In this paper I address the question of whether long-lived gluinos having 
mass less than ~ 1.5 or greater than 3.5 GeV are excluded on other grounds. 
Many experiments which are commonly cited as ruling out gluinos of this 
mass range actually provide only weak limits when one takes account of the 
gluino lifetime. These experiments as well as the most powerful indirect con- 
straints, which are also presently unable to exclude this mass range, will be 
reviewed below. My purpose here is to propose tests which will unambigu- 
ously demonstrate or exclude the existance of light gluinos. 

An inevitable consequence of the existance of a long-lived gluino is the 
existance of neutral hadrons containing them. Generically, hadrons contain- 
ing a single gluino are called -R-hadrons||12|. The lightest of these would 



be the neutral, flavor singlet gg "glueballino" , called R°. There would also 
be i?-mesons, qqg, and R-ba.iyoias,qqqg , with the qq or qqq in a color octet. 
Unlike ordinary baryons which are unable on account of fermi statistics to 
be in a flavor singlet state, there is a neutral flavor-singlet i?-baryon, udsg, 
called S° below. It should be particularly strongly bound by QCD hyperfme 
interactions, and probably is the lightest of the _R-baryons|2], [L3|], even lighter 
than the i?-nucleons. 

The strategy pursued here is to identify production and detection mech- 
anisms for the R° for which reliable rate estimates can be made, so that 
searches which are sufficiently sensitive will definitively rule them out or 
find them. First, we use theoretical arguments to estimate i?-hadron masses 
as a function of gluino mass. Then experiments are proposed to settle the 
question. 

2 i?-hadron mass estimates 



photon energy spectrum, and thus be excluded by CUSB. 



If the gluino is heavier than ~ 3.5 GeV, then the R° and S* will have 
masses approximately equal to the mass of the gluino. For the window 
nig < 1.6 GeV, lattice gauge theory (lgt) should be used to determine the 
hadron spectrum, and hopefully the necessary calculations will be done soon. 
However we can get a rough idea without it, as follows. Let us begin by es- 
timating hadron masses if the gluino is as light as possible. 

If the gluino were massless, the spectrum would be expected to contain 
an unacceptably light Jl4, [15| flavor-singlet goldstone boson associated with 
the spontaneous breaking of the non-anomalous linear combination of quark 
and gluino chiral U(l) symmetries. For three light flavors of quarks the 
non-anomalous axial current isQ 

4 = -m i^^ ~ ffl-r^L - ^A a } . (i) 
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We can obtain a theoretical lower bound on the gluino mass by identifying 
the rf with this pseudogoldstone bosonF[ The flavor singlet pseudoscalar 
which gets its mass from the anomaly would then be identified with a more 
massive state, which will be discussed below. 

If this were the correct description of the 7/, its quark content would be 
reduced by a factor of || ~ 0.7 in comparison to the usual picture. In- 
terestingly, this seems not to be ruled out by existing constraints. Sound 
predictions for the 77', avoiding model dependent assumptions such as the 
relation between F\ and F$, are for ratios of branching fractions to final 
states which couple to the quark component ]16[. These ratios are insensitive 
to the presence of a gluino or gluonic component. Absolute predictions are 
highly sensitive to theoretically incalculable hadronic effects, due to the very 
restricted phase space for the 7/ to decay through strong interactions. This 



The fields appearing in this expression are left-handed Weyl spinors and a sum over 
indices is understood, i labels the three light quark flavors and j and a label the 3 quark 
and 8 gluino color degrees of freedom. 

6 This possibility was suggested in ref. M but not developed in quantitative detail as 
is done here. 



means that rates which could potentially determine whether the rf has a 
30% gluino component, in practice cannot be predicted reliably enough to 
be useful.Q 

Assuming tat the rf is the pseudogoldstone boson connected to the spon- 
taneous breaking of the conserved axial current (|l|), like the K + for J\ = 
~m {^l1^ s l ~ uC l1v sC l p standard current algebra manipulations lead to 
predictions for m^,f^ and m^f^. Taking their ratio, neglecting m u and 
irtd in comparison to m s , and solving for m g leads to: 



Wig ~ 



m s < ss > 
T< AA> 



13 ( m W 



K , 



(2) 



With f v > « Jk this gives m g ~ ll~~^m s . Since the QCD attractive force be- 
tween color octets is greater than that between triplet and antitriplet, < A A > 



is presumably larger than < ss >. Most-attractive-channel argument s [|T7 
suggest that the condensates depend exponentially on the Casimirs of the 
condensing fermions so that since Cs/Cs = 9/4, < AA > could be an order 
of magnitude or more larger than < ss >. Thus pending lattice calculations 
of < AA > or m(r]') as a function of gluino mass and without gluinos, the 
phenomenological analysis should be general enough to include a gluino as 
light as ~ 100 MeV or less. In this case the i?-hadron properties are about 
the same as they would be for a massless gluino. 

If the gluino were massless, the mass of the R° should be 1440 ±375 MeV, 
i.e., about l\ GeV, as follows. Consider supersymmetric SU(3) Yang Mills 
theory. Since supersymmetry in this theory does not break dynamically |T8 



hadrons must fall into degenerate supermultiplets. The massive chiral su- 

permultiplet containing the ++ glueball also contains a h (the lowest gg 

bound state) and two spin-| states, namely the two helicities of the R° (the 

gg bound state). At the classical level this theory has a chiral U(l) phase in- 

variance since the gluinos are massless, like the chiral U(l) of ordinary QCD 

7 A possible way to discriminate is to study the production of the various pseudoscalars 
in J /^> decay. G. Farrar and G. Gabadadze, in preparation. 



with massless quarks. This symmetry is clearly not realized in the hadron 
spectrum, since the R° is degenerate with the massive glueball. Nor is there 
a goldstone boson associated with the breaking of this U(l) symmetry since 
the pseudoscalar gg bound state is also degenerate with the glueball. This is 
not paradoxial for the same reason that in ordinary QCD we can accomodate 
the t]' mass. Namely, the axial U(l) current has an anomaly so that non- 
perturbative effects give the pseudoscalar gg bound state a mass. The chiral 
U(l) symmetry is explicitly broken by quantum effects so that Goldstone's 
theorem is circumvented^. 

Now consider hadron masses in QCD with three light quarks and a mass- 
less gluino. The mass of the ++ glueball is predicted |1| using lattice QCD 



in quenched approximation (i.e., QCD with only gluons and no quarks or 
gluinos) to be 1440 ± 110 MeV. In ordinary lattice QCD, with three light 
quarks but no gluinos, the quenched approximation is commonly taken to 
be valid at the 10 — 15% level for hadron masses^. Since the 1-loop beta 
function for QCD with no light quarks but an octet of gluinos is the same as 
for QCD with three light quarks, one can expect that the error for quenched 
approximation in the supersymmetric Yang Mills theory is also 10 — 15%, 
and that the quenching error with both quarks and gluinos is ~ 15 — 25%. 
If this is so, then a full lattice calculation in QCD with 3 light quarks and 
a massless octet of gluinos would give a mass for the R° of ~ 1440 ± 375 



MeV, where the lattice error of ref. [IS] on the glueball mass was combined 



in quadrature with the estimated error from the quenched approximation, 
taken to be 25% of 1440 MeV to be conservative. As we shall see below, it 
is much more difficult to detect an R° with mass ~ 1065 GeV than one with 
mass ~ 1800 GeV. Thus a lgt calculation which reduced the range of uncer- 



8 It is interesting that supersymmetry relates the mass produced by non-perturbative 
effects through the anomaly to the mass-gap for the glueball coming from confinement, 
suggesting that confinement is essential to the understanding of the mass of the rf even 
in ordinary non-susy QCD. 

9 See ref. |2fl] for a critical discussion of quenched approximation. 



tainty on the mass of the R° would be very helpful, especially if it showed 
we could ignore the region close to 1 GeV. 

In QCD extended by gluinos, the flavor singlet pseudoscaler which gets 
mass from the anomaly is orthogonal to the anomaly- free current (p]) , thus it 
is 70% gg and 30% uu + dd + ss. In the supersymmetric Yang Mills theory 
discussed above, the pseudoscalar gg state which gets mass from the anomaly 
and is degenerate with the ++ glueball would have a mass of 1440 ± 240, 



adding the error in ref. [19| in quadrature with a 15% error for unquench- 



ing. There is evidence for an "extra" flavor singlet pseudoscalar present in 



the meson spectrum in the 1410-1490 region|21, G2l E3j, which has a large 



coupling to gluons|TlJ- If confirmed, it is an excellent candidate to be the 
pseudoscalar whose mass comes from the anomaly, in the very light gluino 
scenario. 

To recapitulate, we have seen above that from purely theoretical con- 
siderations we can at present only rule out R° and S° masses below about 
1100 MeV. Having the lightest possible masses requires both the r{ and extra 
pseudoscalar meson in the 1410-1490 MeV region to have large gluino com- 
ponents, but increasing the gluino mass to ~ 700 MeV allows one to return 
to the conventional phenomenology for the rf and interpret the extra state 
as a simple r]g, the lowest-lying gg bound state. If gluinos are much heavier 
than this, one needs another explanation for the extra state in the 1410-1490 
region. 

3 Existing Experimental Limits 

From the CUSB experiment, we inferQ that the r) g does not lie in the 3-7 
GeV range, so that the gluino would not be in the ~ 1.5 — 3.5 GeV range. 
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See footnote on the first page of the present paper. 



In order to compare to limits from other experiments searching for R 0, s, we 
shall convert this limit to an effective gluino mass using the relation 

m(R°) = 0.72(1 + e-~) + m- g {l - e~ m §), (3) 

with all masses in GeV. This is actually just a convention for making the 
figure, but is physically reasonable in that it yields the nig = result of the 
previous section and in analogy with mesons made of one light and one heavy 
quark associates an additive confinement energy of about half the mass of a 
light-quark-meson (here, of the ++ glueball whose mass is ~ 1.44 GeV) to 
the light constituent (here, the gluon) of a light-heavy composite. 



In another quarkonium decay experiment, the ARGUS group [24| looked 
for events in which T' — > 7 + Xb(l 3 -Pi), followed by x&(l 3 Pi) — > ggg, with 
one of the final i?-hadrons decaying in a distance of 1-60 cm from the e + e~ 
interaction point. From the absence of such events at a level predicted by 
pQCD they concluded that gluinos in the mass range 1-4.5 GeV do not exist 
in the lifetime range to which they were sensitive. However perturbative QCD 
overestimates the branching fraction Xb{^ 3 Pi) -^ 999 for very light gluinos, 
since it fails to include the effect of the substantial reduction in phase space 
arising from the minimum invariant mass of a pair of R 0, s being about 3 GeV, 
even when the gluino is massless (see section ||). To determine whether the 
experimental sensitivity extends to a gluino mass as low as 1 GeV as stated 
in ref. p4j , the experiment should be reanalyzed using a more realistic model 
of the branching fraction for Xb{^ 3 Pi) ~^ 999 m the non-perturbative portion 



of phase space. The ARGUS results, taken from Fig. 4a of ref. p4 



plotted on the figure using the above function to convert from their quoted 
gluino masses to a common R° mass. For the largest masses no conversion 
is used, in order not to make the non-sensical claim that they can exclude 
R 0, s which cannot be kinematically produced. 

The best constraints beyond CUSB and ARGUS for long-lived gluinos in 
the radiatively-generated range of up to O(30) GeV come from searches for 
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new neutral particles. Gustafson et al.|25|| searched for new hadrons with 
lifetimes greater than 10~ 7 sec, using time-of-flight in a 590m long neutral 
beam at FNAL. On account of timing and energy resolution limitations, they 
were capable of distinguishing a particle from a neutron only if its mass was 
greater than 2 GeV. From the limits of Gustafson et al, Dawson, Eichten and 
Quigg|| (DEQ) concluded that gluino masses in the 2-4 GeV range could be 
excluded. This experiment is therefore consistent with CUSB and Bernstein 
et al (see below), and for r g > 10~ 7 sec extends the lower end of the excluded 
mass range to 2 GeV as shown in the figure. 

The experiment of Bernstein et al.[^(| places an upper bound on the 
production cross-section of a neutral hadron produced in 400 GeV proton 
collisions, with mass in the range 1.5 — 7.5 GeV, which decays with a lifetime 
(10~ 8 — 2 x 10" 6 ) sec to a 2- or 3-body final state containing a charged 
hadron. They find E^-\ 90 o<5 x 10~ 35 (Gc ^ /c 3) f° r mass of 1.5 GeV, and 
<3 x 10~ 32 ( Gc '^2/ c 3) for 7.5 GeV, taking the most sensitive lifetime value 
of 3 x 10~ 8 sec. Typical decays would be R° —>■ lsp + vr('s) and S° — > 
lsp + A + vr('s) or S° — > lsp + N + K + vr('s). Since the S° has baryon 
number +1, it would be expected to be produced mainly in the forward 
direction rather than at 90° where the experiment was done, so is not directly 
constrained by this experiment. However this experiment does constrain the 
possibility of _R 's. For the light end of the mass range, a reasonably good 
analog process which should be even more OZI-suppressed, is pp — > pX 
whose invariant cross section is ~ lO -27 ,^™ ^ [p7|| , for similar kinematics. 



For a gluino mass of 3.5 GeV or larger, it is legitimate to use perturbative 
QCD (pQCD) to compute the expected rate, as a function of gluino mass. 
This was done in tree approximation by DEQ|| for m g = 3 GeV. They 
predicted an invariant cross section of ~ 10~ 28 tq^tt^] for p± = 0. We can 
very crudely estimate the cross section for production of a gluino of higher 
mass but p± = by noting that the cross section is mainly dependent on 
the combination m 2 + p\. The DEQ prediction for m = 3 GeV and p± = 4 



GeV is ~ 3 x 1Q~ tqw^ (see Fig. 44), which is the same as the Bernstein 
et al limit for m = 5 GeV and p± = 0. Thus the Bernstein et al limit very 
roughly rules out R 0, s with mass less than 5 GeV. The range of lifetime 
sensitivity corresponding to the cross section limits of Fig. 4a is shown in 
Fig. 4b, for m = 3 GeV where it is ~ 2 x 10~ 8 — 2 x 10~ 7 sec. Since 
for a fixed production rate the detector sensitivity depends mainly on 7/3r, 
and 7 ~ mT 1 , the range of maximal sensitivity will shift upward, roughly in 
proportion to m, for m > 3 GeV. The range excluded by Bernstein et al is 
shown in the figure. It is the upper elongated region ending at 5 GeV. 

The limits could in principle be somewhat tightened if there are charged 
i?-hadrons which decay only weakly to the R° or S°, e.g., Rk+ — » i?° + 7r + . 
This will be the case if the mass gap between charged .R-pions and i?-kaons 
and the R°, and between charged i?-baryons and the S°, is greater than 
the mass of the corresponding kaon or pion. Lattice calculations of the R- 
hadron mass splittings as a function of gluino mass are badly needed here. 
Bag model predictions for .R-hadrons cannot be trusted since parameters 
fixed to fit the ordinary hadrons may not be applicable to .R-hadrons, and 
furthermore bag model estimates have not been been reliable for the glueball 
spectrum. Nonetheless old bag model estimates |S8|, |2], [T3| suggest that for 



some parameters there may not be enough phase space for Rk — > K + 
R° or Rn — > K + S°. Thus a search for charged i?-hadrons is worthwhile 
even though a null result would not exclude gluinos. Note that there is no 
relation between the lifetimes of the R° and S° and lifetimes of charged R- 
hadrons, since the latter decay to the R° and S° through conventional weak 
interactions and would be expected to have a lifetime comparable to weakly 
decaying hadrons of a similar mass, i.e., 10~ 10 — 10~ 13 sec for masses in the 
range 1-5 GeV. Briefly, the experimental constraints would be: 



Cutts et al|29| use time of flight to exclude lifetimes greater than ~ 
(2 — 5) x 1CT 8 sec, for charged particles with masses in the 4-10 GeV 



range. 



Bourquin et al[^] search for decaying particles in the CERN hyperon 
beam, extending the excluded range for new charged particles to cover 
the 2-4 GeV mass range, for lifetimes of order 1CT 9 — 10~ 8 sec. 

Charged R— hadrons having mass of the same order of magnitude as 
the D or B mesons must have a lifetime too short or long to decay in 
vertex detectors used to measure D and B lifetimes. 

There is a CDF limit on the existance of charged hadrons having 
7T>1CT 7 sec PHI, but it only addresses masses greater than 50 GeV 



because the present detector has time resolution at the nanosecond 
level. 

Otherwise the constraints on charged i?-hadrons are poor and the cover- 
age is surprisingly spotty. It must be reemphasized, however, that even if 
strong-interaction-stable i?-hadrons exit, one cannot immediately apply these 
experimental constraints on the allowed regions for their mass and lifetime 
to the limits in the figure, because there is no direct relation between the 
lifetime of a charged .R-hadron and that of the R°. 

If the gluino lifetime is long because the squark mass is much larger than 



m\y, then beam dump experiments! 12, B3, p3, p4l B5f, which look for the 



reinteraction of the lsp in a neutrino detector, become ineffective because 
the lsp cross section falls as ( ^p- J . Even if the lsp cross section is not too 
small, the gluino must decay before losing energy in the dump, e.g., in 10 



cm in the Ball et al FNAL beam dump experiment [p3| , p6| , i.e., requiring 
a lifetime ^5 -^y lO" 11 sec. Likewise, the BEBC experiment ||34|| observes 
that if Tg > 5 x 10~ n sec the gluino decay does not occur before interaction, 
"severely degrading the photino flux reaching our detector." For massless 
photino they model this effect, but in general beam dump experiments need 
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to be analyzed in terms of the three parameters m§, a\ sp and r g . Beam dump 
experiments cannot be used to exclude regions of the gluino mass-lifetime 
plane without further assumptions which are not in general appropriate to our 
case, except for gluinos with lifetimes shorter than about 5 x 1CT 11 sec. The 
HELIOS experiment J35[1 explicitly addresses direct production of WINP's, 



and not long-lived gluinos, since it requires that no energy degradation occur 
in the dump. 

The possibility of large gluino mass is at present only addressed by collider 
missing energy searches that detect the existance of a gluino which decays 
inside the apparatus with a substantial portion of its energy going to the lsp 
which is very weakly interacting and escapes. Indeed, this is the classic gluino 
signal J37| . The CDF missing energy searchjjlj is sensitive to gluinos which 



decay within about 1 meter of their production, i.e., having — -r g <3x 10 -9 
sec. They require the missing transverse energy to be greater than 40 GeV. To 
get a very rough idea of their regime of sensitivity (which could be determined 
more accurately by modeling the energy spectrum of the produced gluinos) 
we can take as a typical event the case in which the gluino is emitted at 45° 
and assume it decays giving 1/3 its energy and momentum to the lsp which 
escapes with the minimal transverse energy to satisfy their cuts. In this case, 
the actual energy of the decaying gluino would be 3 ■ V2 -40 = 170 GeV, 
ignoring gluino, quark and lsp masses. Thus gluinos with lifetimes longer 
than about 2 x 10~ n (y^v) sec ' wou ld not be efficiently detected in the 
CDF search. They do not investigate masses lower than 20 GeV, where they 
lose efficiency on account of the acoplanarity and missing E? cuts. 

The UAl missing energy search |38|] claims to be sensitive enough to ex- 
clude masses as low as 4 GeV. Although a gluino lifetime is not included 
in their efficiency monte carlo, they state that they believe they are fully 
sensitive to gluinos whose lifetime is shorter than 10~ 10 sec. This agrees with 
the crude estimate given above for the CDF experiment, for rrig = 5 GeV. 
Nonetheless, especially for the high mass end of the UAl experiment, a monte 
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carlo is needed to know the lifetime sensitivity as a function of mass. For 
simplicity, and to be conservative, we will use the estimate 2 x 10~ n ( x q g v ) 
sec for both CDF and UA1. 

To summarize this section, gluinos in the mass range ~ 1.5 — 3.5 GeV are 
absolutely excluded (CUSB). Lighter gluinos are allowed, as long as the R° 
lifetime is not in the range 2 x 1CT 6 — 10~ 8 sec if the R° mass is greater than 
1.5 GeV (Bernstein et al), or the range > 10~ 7 sec if its mass is greater than 
2 GeV (Gustafson et al). Gluinos with mass around 4 GeV or above, must 
have a lifetime longer than about ~2x 10 -11 (y^fy) sec (UA1,CDF), with 
the ranges > 1CT 7 sec (Gustafson), 2 x 1CT 6 — 10~ 8 sec (Bernstein et al) and 
~ 1CT 10 sec (ARGUS) ruled out for masses in the vicinity of 4-5 GeV. The 
figure is an attempt to summarize these results, combining experiments which 
report results directly in terms of m(R°) with those characterized by limits 
on nig by use of eqn. (|3]). Given the primitive nature of eqn. (|3|) and the 
±375 MeV uncertainty on the R° mass when the gluino is massless (section 
|2|), as well as the very rough methods used to extract the ranges of mass and 
lifetime sensitivity for the various experiments, a > 20% uncertainty should 
be attached to all the boundaries shown in this figure. 

4 Theoretical Comments: Gluino Lifetime and 
Production Estimates 



How natural is it from a theoretical point of view for an R° in the mass 
range 1.5 — 2.5 GeV to have a lifetime longer than 2 x 10~ 6 sec, or for an R° 
with mass >5 GeV to have a lifetime longer than ~ 2 x 10 -11 ( 1 Q ( fv ) sec ^ 
For the higher mass range the R° and gluino lifetimes can be taken to be 
approximately the same, since for a relatively massive state one can ignore 
the effects of confinement on the overlap of the initial and final states, and 
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the modifications to phase space from the hadron masses. For the low end of 
the range, if the lsp mass is low compared to the gluino mass, one could either 
argue by analogy to known hadron decays |12[] or, following Franco |39|, take 



the R° lifetime to be that of a gluino of ~ | of its mass. For the interesting 
case that the lsp mass is a significant fraction of m(R°), tools have not yet 
been developed which allow us to reliably estimate the resultant suppression 
in the decay rate. 

The decay rate for an unconfined gluino to decay to the lsp and a qq 
pair can be obtained as follows. In general, the lsp is a superposition of the 
fermionic partners of the neutral SU(3) and U(l) gauge (w3-ino and bino) 
and Higgs bosons (higgsinos). However it is shown in ref. that when 
gaugino masses are all radiatively generated the higgsino component of the 
lsp is in fact less than 1% in amplitude. Thus we can approximate the lsp 
wavefunction as cos9\b > +sin9\ws >. The decay rate of the gluino assuming 
the lsp to be a photino was given in ref. [fJ0| , so we need only replace e 2 
appearing in their expression by ( J™f )[I Z + z^} 2 , where z = *°"^ , and 
average over left and right handed contributions. Thus the total rate for 
gluino to decay to the lsp and a uu, dd, or ss pair, ignoring the quark 
masses is: 

_ q s q em (l-fz+z 2 )m| sinQ „ 
1 9 1287rM s 4 g v sinO w > W 

[(1 - y 2 )(l + 2y- 7y 2 + 20y 3 - 7y 4 + 2y 5 + y 6 ) + 2% 3 (1 -y + y 2 )log(y)], 

where y = ^ i££ . We have taken M S jJ = M^ = M sq for simplicity. The 
9 dependent factor ranges from 1, for a light neutralino in the low-// region 
where 9 ~ 9\y, to (cos9w)~ 2 for a heavy neutralino in the high-/z region where 
cos9 ~ 1. Thus for a rough estimate we take this factor to be 1. We also 
take a s ~ 0.1 and a em = 1/128, since the relevant scale is the squark mass. 
Then, for instance with a massless lsp, the squark mass must be greater than 
~ 2 TeV for a gluino with effective mass in the 1-1.5 GeV mass range to have 
Tg > 2 1CT 6 sec. If instead the lsp mass is 90% of the gluino effective mass, 
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the squark mass must only be greater than about 200 GeV. For a gluino of 
mass 5 GeV, the UAl bound is most relevant. For lsp mass of zero or 0.9rrig 
one finds that the squark mass must be greater than 1 TeV or ~ 130 GeV, 
respectively. These squark masses increase to 6 TeV or 670 GeV for a 15 
GeV gluino. As shown in ref. ||, when gaugino masses arise radiatively, 
these conditions are naturally accomodated in much of parameter space. 

It is also worth noting that absolute stability is a real possibility for the 
S°, since the mass difference between it and the lsp must be greater than 
the sum of proton and electron masses for it to decay. If it binds to nuclei, 
this would be ruled out experimentally by the sensitive searches for exotic 
isotopes, at least for some mass regions||l|. However one would expect a 



repulsive, not attractive, interaction between a nucleus and the flavor-singlet 
R° or S° , since the intermediate state created when they exchange mesons 
with a nucleon has a much higher energy^. 

Anomalous signals in extensive air showers and underground muons seem- 
ingly coming from Cygnus X-3 are consistent with the intermediate particle 
being a neutron, except that the neutron decays too quickly to make the long 
tripQ. Long-lived R 0, s were investigated [^], but discarded^] on account 



of the mistaken belief that they would imply a long lived charged i?-proton 
which is ruled out by, e.g., ref. J4l| . If the present quiet of Cygnus X-3 is 



only a cyclical phenomenon and such events are observed again in the future, 
an S° interpretation should be seriously considered. 

Turning now to cross section calculations, I am not aware of any re- 
cent pQCD calculations of gluino production at a hadron collider, except for 
very massive gluinos. The old analyses P, |45| should be updated, making 
an attempt to estimate the uncertainty in the gluino distribution, as well 



11 Unlike the binding of a nucleus where exchange of mesons between pairs of nucle- 
ons, each of which can absorb or emit an I = 1 meson and remain a nucleon, leads to 
intermediate states close in energy to the original state. 

12 See, e.g., ref. |42| for a summary. 



14 



as including 1-loop corrections which have proved very important for ordi- 
nary pQCD predictions. From deep inelastic and Drell-Yan experiments, the 
quark and antiquark distributions are reasonably well fixed. Direct photon 
production gives information on the gluon distribution function, so the mo- 
mentum sum rule then provides some constraint on the gluino distribution. 
The naive argument p6| which leads to behavior (1 — x) 7 for the sea-quark 
distribution functions at large-x, leads to the same behavior for the gluino 
distribution function. Since the R°, rj', and rjg masses are so much larger than 
pion masses, one would expect that the low-Q 2 gluino distribution functions 
are smaller than those of the sea-quarks. However since the 1-loop beta- 
function for gluinos is the same as for 3 flavors of light quarks, the gluino 
distribution function evolves as rapidly as all three quarks together, so a light 
gluino would become an important component of the nucleon at larger Q 2 . 

Although the gluon and gluino distribution functions are individually dif- 
ficult to determine well, without assumptions as to their functional form for 
the entire x range, their sum is much better determined^. Since both gluons 
and gluinos give rise to gluino jets, the actual prediction for R° production 
is relatively stable. If the existance of gluinos were established, the ratio of 
events with 1 and 2 R 0, s would allow the ratio of gluino and gluon distribu- 
tions to be constrained. Demanding consistency of pQCD predictions with 
observed jet production may also allow the gluino distribution function to be 
further constrained, since the amplitudes for gluinos to produce jets differs 
from those for quarks or gluons to produce jets. 

5 Indirect Evidence Regarding Light Gluinos 

For years it has been recognized that in principle the running of a s is 



13 Comparably to the determination of the gluon distribution function, when the gluino 
possibility is ignored. 
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sensitive to the presence of gluinos0. In deep inelastic scattering experi- 
ments the ambiguity introduced by higher twist contributions is too large 
to allow one to decide between QCD with and without gluinos0. Gluinos 
modify the e + e~ annihilation cross sections only in order ot 2 s , by providing 
an additional source of 4-jet events and making virtual corrections to 2-jet 
events. The possibility of infering or excluding gluinos directly from LEP 
event characteristics was discussed in ref. ||49|| , where the sensitivity to the 



as-yet-uncalculated 1-loop corrections was shown to be too great to allow one 
to decide between ordinary QCD and QCD with massless gluinos^. 

The reason that it is generally difficult to discriminate between QCD with 
and without gluinos is because adding gluinos to the theory modifies it in 
competing ways which tend to cancel. For instance the value of a s at LEP is 
obtained by fitting QCD predictions for various aspects of event shapes and 
extracting the value of a s which gives the best fit . Gluinos are an additional 
source of 4-jet events, but at the same time a s runs more slowly when there 
are gluinos. This means that, for a given value of a s (Mz), the typical value 
of a s {Q e ff) in multi-jet events is lower than it would be for QCD without 
gluinos, which tends to reduce the number of multi-jet eventsQ. 

Just as the effects of gluinos tend to cancel at LEP, one cannot simply 



14 For the first discussion of this, see ref. [E7[ 

15 Thus comparison of the values of a s from deep-inelastic scattering and Z° decay are 
inconclusive, although suggestive Q. 

16 More recent articles on this subject have come to the same conclusion |p0 51 1. 



17 Ref. |32| correctly emphasized the need to extract a s from data with and without 
gluinos before evaluating the consistency of the running between different energy scales, 
with and without gluinos. However that analysis only includes the virtual corrections to 
the running of a s in LEP events and not the effect of real gluino jet production which 
is of the same order, so is incomplete. Bryan Webber and I (unpublished) tried to see 
if we could find some systematic preference of the LEP data for QCD with and without 
gluinos by looking at the entire menagerie of quantities from which a s is extracted. We 
found that the only region in which there was a significant difference in predictions with 
and without gluinos is precisely the region in which hadronization is most important, and 
for which the 1-loop corrections to the 4-jet cross section (which are not yet available) are 
crucial. 
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say that the number of jets predicted at the Tevatron will be increased by 
such-and-such an amount, since if there are light gluinos they will be present 
in the hadron structure functions and will use some of the "room" for gluons, 
reducing the production of conventional jets to some extent^. To address 
the possibility of light gluinos by their effects on jets or the running of a s , 
one must a) compare predictions for actual experimental observables with 
and without gluinos and not try to compare derived quantities such as a s 
and b) fully incorporate gluinos into the analysis, including their effects on 
distribution functions. 

Recently, there have been a number of attempts to make the kind of care- 
ful analysis which would be necessary to obtain reliable indirect information 



of the possibilty of light gluinos. Ref. [53] used theoretical predictions for the 



hadronic branching fractions R T and Rz with and without gluinos, to extract 
a s (m T ) and a s (mz) with and without gluinos, then checked whether the run- 
ning of a s between these values was consistent with what QCD predicts, with 
and without gluinos. The main difficulty with this approach is the issue of 
how to treat the effect of a light or massless gluino on r decay, and also the 
question of the validity of neglecting non-perturbative contributions of order 
■£p as is done in ref. [|54| . The latter issue is discussed in ref. | j55|| , where it 



is argued that unless the validity of neglecting ^ corrections is established, 
the error should be taken to be twice that assigned in the "nominal" case of 



ref. |p3| , i.e., that Rth = 2 is appropriate for the ref. |k| analysis. With 
respect to the former issue, since the invariant mass of an R° pair and of 
the rjg is too large to contribute significantly to r decay, independent of the 
gluino mass, the gluino contribution should be neglected when determining 
a s (m T ), as is done for the charm quark. This can be implemented^ by using 

18 See ref. |49] for a discussion of the difficulty of using hadron collider jet cross-section 
characteristics to infer or exclude the existance of long-lived gluinos given the present level 
of theoretical precision. 

19 M. Schmelling and R. St. Denis, private communication. 
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an "effective" gluino mass >m T /2 when using their fig. 2. One then finds 
that even at 90% confidence level there is no excluded region of gluino mass 
from this analysis when Rth = 2. 

For other attempts to study indirect evidence for light gluinos see, e.g., 
refs. [^6], [57 1 . 



6 Proposals for Experiments 



Now let us turn to the question of how to establish or rule out the exis- 
tance of new light hadrons, R° or 5°. One method, proposed years ago[Q], 
is to look for exclusive reactions such as K~p — > R°S°, followed by elastic 
scattering of the R° and S° off protons. With accurate measurements of the 
R° and S° production angles, and measurement of the recoil proton momenta 
in the secondary R°p — > R°p and S°p — > S°p scatterings, there is in principle 
one more equation than unknowns and the masses of the R° and S° can both 
be determined. Using a hydrogen bubble chamber would seem to work nicely 
for observing the initial and secondary scatterings, but a high efficiency for 
identifying if£'s and neutrons would be desirable to reduce background, so 
this may not be the optimal approach. The interaction lengths of the R° and 
S° are probably somewhat shorter than for ordinary mesons and baryons, 
on account of the greater color charge of the gluinos as compared to quarks 
and on account of the S° having 4 constituents rather than 3 for a normal 
baryon. The candidate events should show a threshold behavior consistent 
with the measured R° and S° masses, which would corroborate the validity of 
the overall picture. Note that this experiment is sensitive to gluinos with any 
lifetime long enough that the R° and S° rescatter before decaying, so that it 
is complementary to the experiment of Bernstein et al. and sensitive to lower 
masses than Gustafson et al. However this method has two important weak- 
nesses: First, the cross section may be very small, since one is asking for a 



very exotic final state to be produced in an exclusive mode. Second, it is not 
possible to reliably calculate the cross section so that one cannot establish a 
level of sensitivity adequate to definitively exclude the phenomenon. Unfor- 
tunately it is also a demanding, single-purpose experiment and theoretical 
prejudice has favored heavy gluinos, so that experimenters have not looked 
just to see if something might be there0. 

Here I propose other experiments which also do not rely on observing the 
decay of the R° or S° and are thus able to rule out or observe long-lived 
gluinos, but which do not have the difficulties of the one discussed above. 
Except in the forward direction, we expect that S° production is much smaller 
than R° production, so let us ignore S°'s for simplicity. By working at high 
energy, exotics can be produced relatively easily and inclusive cross-sections 
can be reliably computed from perturbative QCD in appropriate kinematical 
regions. The cross section for producing R 0, s is essentially just the gluino 
jet cross section, since all gluino jets end in an R° (or, rarely S°) because 
other i?-hadrons eventually decay to theseQ. The gluino-jet cross-section is 
approximately 10%||49| of the total jet cross section, so that it is actually 



quite common for a Tevatron collider or fixed target event to contain an 
R° pair. p t cuts can be imposed to insure that perturbative QCD event 
generators can be reliably used to compute the expected rate, even for light 

20 Recently, Carlson and SherBS] proposed searching for the decays of gluinos following 
their photoproduction at CEBAF. This is an excellent experiment, since something may 
be found. However it does not satisfy the present criterion of being useful for excluding a 
light gluino, since the relatively low invariant mass range which can be probed at CEBAF 
means that the non-perturbative effects of R° and rjg masses will suppress the signal 
and the calculations of the production rates are therefore not sufficiently reliable to allow 
exclusion. They report results for the effective gluino mass being taken to be 1 and 1.5 
GeV, and the dramatic rate of decrease with effective gluino mass reflects the sensitivity 
to this effect. To obtain reliable inclusive cross sections for production of light particles 
from pQCD, one must impose a p™ m cut. The event rates they quote are so large that 
this may be possible, but as long as their signal is the decay of the gluino, the proposed 
experiment can only be used to examine a limited lifetime region. 

21 Or very rarely, gluinos from independent jets can annihilate, but at this order one 
must also consider jet evolution which produces gluinos. 
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gluinos. Showing that there are no such events at a level of 4a below the 
prediction, would then convincingly rule out the existance of these gluinos^. 

Basically the idea is an outgrowth of the suggestion of ref . , but sacrific- 
ing the additional constraints of exclusive production in favor of the higher 
rate and reliable calculability of high energy inclusive production. A high 
energy beam from an accelerator is incident on the primary target. This 
produces a neutral beam containing neutrons, kaons, hyperons, and possibly 
i?°'s and S 0, s. This beam illuminates a secondary target in which an elastic 
scattering R°p — > R°p may occur. Measuring the momentum of the recoil 
proton and the angle of the produced R° (by observing its interaction, which 
need not be elastic) gives enough constraints to solve for rriR, if indeed the 
reaction is elastic. Knowing the visible energy of the final particles in the 
secondary scattering of the produced R° can help choose between multiple 
solutions and help discard events in which the primary scattering is not elas- 
tic.Pl Of course the background due to other reactions, especially n p — > n p 
or K\p — > Kip, or inelastic scattering, will be quite severe even after vetoing 
on extra charged particles and 7r°'s, so excellent resolution is crucial. 

Timing could be used to measure p/E of the incident neutral. With this 
information, one would have an over-constrained system of equations without 
relying on the secondary scattering being elastic and one could verify that 
the initial reaction was indeed R°p — > R°p as well as determine the R° mass. 
If the R° is sufficiently heavy, one can get adequate resolution with nsec 
accuracy using the beam buckets without being forced to put the secondary 
target so far away that the loss of solid angle would be intolerable^]. Modern 

22 Care must be taken to realistically estimate the theoretical uncertainty, including that 
from the distribution functions and neglect of higher order corrections to the partonic 
scattering amplitudes, which in ordinary QCD have proven to be larger than originally 
estimated. 

23 I am grateful to T. Devlin for making these points. 

24 Keeping the distance between the two targets as small as possible is also desirable 
from the standpoint of being sensitive to relatively shortlived i?°'s as well. 
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0(10) psec timing could allow the lower mass regions to be investigated, 
except that it requires tagging the initial R° production event, so entails a 
reduction in rate. Detailed monte carlo simulation is needed to determine 
whether it is possible to cover the very-light gluino regime, where the R° may 
be difficult to distinguish from a neutron. With many events, a discrepancy 
between the observed and expected event characteristics such as angular 
distribution and rates would be a useful diagnostic. Another handle for 
some range of R° lifetimes would be a distance dependence of the anomalous 
events. 

In the above discussion I focussed on the process R°p —>■ R°p for iden- 
tifying the R°. It is the most attractive option from a theoretical point of 
view since its cross section is easiest to estimate^]. If one's goal is to try 
to unambiguously exclude light gluinos, then one must use reactions which 
can be estimated with some confidence both to produce and to detect them. 
However if one wants the most effective way to discover light gluinos if they 
exist, one can consider other detection reactions such as R°p — > R p i]' or 
R°p — > K + S°, whose signature may be much more distinctive! 26 ]. In the res- 
onance region, such cross sections can be very large. Further work is needed 
to try to estimate them. 

A setup such as KTeV, where the distance between primary and secondary 
targets (the regenerator) is 120m and the typical energy of the long-lived 
neutrals is about 100 GeV, would be mainly sensitive to lifetimes longer 
than ~ 4 10~ 9 sec. Thus if it can be used for this purpose, it will be able to 



25 The optical theorem relates the forward elastic cross section to the total cross section. 
Above the resonance region one would expect a(R°p) ~ a(irp) ~ cr(pp) since the confine- 
ment scale rather than the color charge of the valence constituents, seems most important 
in determining the size of a system of light, relativistic quarks or gluons or gluinos. Using 
lattice gauge theory, it might be possible to measure the color charge radius of the R° 
or at least its ratio to that of the pion or nucleon, to improve upon this crudest possible 
estimate. Or one could use information from lgt on glueball masses to try to constrain a 
bag model for color octet constituents, and then determine their radius. 

26 I am indebted to W. Willis for emphasizing this point. 
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probe a large part of the interesting lifetime range. 

In a collider experiment, pair produced heavy gluinos would radiate glu- 
ons and light quarks to produce jets containing ordinary hadrons and an R°. 
For sufficietly heavy R° and good timing capabilities, one could in principle 
detect the time delay p/E for the late-arriving neutral particles to deposit en- 
ergy in the calorimeter. Assuming each of them to be an R which stopped in 
the calorimeter, producing very light particles, the energy it deposited in the 
calorimeter would be roughly of the same magnitude as p of the R°. Knowing 
p and p/E, one could solve for mj. A detailed study of the conversion of 
an R 0, s momentum to the energy deposition in the calorimeter (in particular 
the extent of the fluctuations to be expected), is needed to see if this method 
is feasible in practice. Another way that the production of a pair of heavy 
long-lived gluinos might be infered in principle, would be to search for events 
in collider experiments in which fitting energy and momentum conservation 
at the jet level requires two of the jets to be given a large mass. 

7 Summary 

As is shown in ref . || , if gaugino masses are generated by loop effects, the 
gluino and lsp masses will be in the range from ~ 100 MeV to < 30 GeV if 
the SUSY and ew symmetry breaking scales are < 10 TeV. Furthermore, in a 
substantial part of parameter space the lsp is near in mass or heavier than the 
gluino, so that long gluino lifetimes are natural. The phenomenology of such 
light, long-lived gluinos is the subject of the present paper. Some aspects of 
the phenomenology of the associated lsp are discussed in ref. . A very light 
gluino (mass of order a few hundred MeV or less) is particularly attractive 
since it emerges naturally when dimension-3 SUSY breaking operators are 
absent from the low-energy theory, as is the case in hidden sector dynamical 



SUSY breaking with no gauge singlet s|o!|. Consideration of the pseudoscalar 
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spectrum is shown to imply that the gluino mass must be greater than ~ 
100 MeV. A very light gluino would lead to new hadrons, the R° (gg) and 
S° (udsg), with masses around \\ GeV. Experiments to definitively rule 
out or discover them are possible but very challenging. Existing direct and 
indirect experimental constraints are reviewed and found not to address the 
most interesting scenarios. Experiments directed at the higher mass range 
are also mentioned. 
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Figure 1: Experimentally excluded regions of m(R°) and T~ g . Horizontal axis 
is m(R°) in GeV beginning at 1.5 GeV; vertical axis is Loqiq of the lifetime 
in sec. A massless gluino would lead to m(R°) ~ 1.4 ± .4 GeV. ARGUS and 
Bernstein et al give the lightest and next-to-lightest regions (lower and upper 
elongated shapes), respectively. CUSB gives the next-to-darkest block; its 
excluded region extends over all lifetimes. Gustafson et al gives the smaller 
(mid-darkness) block in the upper portion of the figure; it extends to infinite 
lifetime. UAl gives the darkest block in the lower right corner; it extends to 
higher masses and shorter lifetimes not shown on the figure. 
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